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r~| , Nuclear suppression of charmoniuni production in proton-nucleus interactions is 

poorly understood, what restrains our attempts to single out unusual effects in 
heavy ion collisions. We develop a phenomenological approach, based on the light- 
cone dynamics of charmonium production, which has much in common with deep- 
ed I inelastic scattering and Drell-Yan lepton pair production. The key observation is 
the existence of a soft mechanism of heavy flavour production, which scales in the 
quark mass and dominates shadowing corrections and diffraction. It naturally ex- 
plains the surprisingly strong nuclear suppression of J/^ at large Feynman-Xi? . The 
low-Xi? region is subject to a complicated interplay of hard and soft mechanisms. 
With evaluated parameters we nicely describe available data on charmonium pro- 
duction in proton-nucleus collisions. Using these results we predict a new process, 
diffr active production of charmonium on a nucleon target, which fraction in the 
total production rate of charmonium is evaluated at 12 %. 
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1 Introduction 

Production of charmonium in heavy ion collisions is believed to be a sensitive probe for 
new physics. In practice, however, hand- waving motivations are not sufficient, and one 
desperately needs a reliable base line. Unfortunately, present understanding of dynamics 
of hadroproduction of charmonium off protons and nuclei is far from satisfactory. Nuclear 
effects are poorly understood, particularly, there is still no reasonable quantitative expla- 
nation of surprisingly strong nuclear suppression of J/\E' production at high Feynman xp- 
It is too adventurous to predict or to interpret the observed nuclear effects in heavy-ion 
collisions, having no idea about what is going on in proton-nucleus collisions. 

I do not pretend in this short talk to establish a reliable baseline, but only want to 
make another step towards it. 

This talk is based mostly on the recent unpublished results obtained in collaboration 
with Jorg Hiifner, as well as on our previous publications, cited in the text. 

2 Poor man's formula or how long does it take to produce a cc 
pair 

There is a wide spread believe, adopted in most recent analysis of data on J/^ production, 
that the cc pair is produced momentarily at the point of interaction of the projectile hadron 
with a bound nucleon, and then it propagates through the nucleus and eventually escapes 
it. This point of view is inspired by the usual perturbative treatment of cc production 
in the rest frame of the charmonium. However, in the rest frame of the nucleus the 
production time is subject to Lorentz time dilation (or the nuclear thickness is contracted 
in the charmonium rest frame). When the production time is longer than the mean 
internucleon distance in nuclei (~ 2 fm) one cannot say anymore that the cc pair is 
produced on one concrete nucleon. This can also be understood as coherence between 
the waves produced at different nucleons. For this reason the production time is usually 
called coherence time (length). 



The time scale of charmonium production has two extremes. In the low energy limit 
the coherent time tc <^ d, where d is the mean internucleon distance, one can use a simple 
(usually called a conventional approach) formula for nuclear suppression, 
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Here the mean nuclear thickness (T) = A^^ J d'^b T'^{b), where T{b) = Tz{b)\z^^oo and 
Tz{b) = J^ dz p{b, z). The nuclear density p{b, z) depends on impact parameter b and lon- 
gitudinal coordinate z. The latter approximation in (|l|) uses smallness of the charmonium 
inelastic interaction cross section, afj^ T{b) <^ 1. 

Although a cc pair, rather than the \E', is produced and propagates through the nucleus, 
in the low energy limit one may think that the charmonium wave function is formed 
instantaneously. We come later back to this point and formulate an exact condition for 
that. 

This simple formula corresponding to the low-energy limit is nowadays widely used as 
a phenomenological basis for nuclear effects in charmonium production. The absorption 
cross section afj^ is usually treated as an unknown parameter. 

In the high energy limit, tc ^ Ra, a fluctuation, containing the cc pair lives much 
longer than the nuclear size. The nuclear suppression factor has in this case a form similar 
to that for quasielastic scattering |jl|], 

ShA\L»RA = -. / d^^ T{b) exp[-afree T{b)] ^ 1 - ajree {T) . (2) 
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I would like to draw attention to the fact that in the approximation af^ T{b) ^ 1 the 
shadowing term in (|I]) has a factor 1/2 compared to that in (Q). This is an explicit 
manifestation of the space-time pattern of interaction: at high energy the mean length of 
path of the cc pair in the nucleus is twice as long as at low energy. 

Let us postpone to the next section interpretation of the freeing cross section a free, 
but for now we concentrate on the problem of the coherence length. The lifetime of a 



fluctuation of tlie projectile hadron containing the cc pair is given by the energy denomi- 
nator, 

_ 2Eh 
' " Mj, - ml ' ^'^ 

where Eh and rrih are the energy and mass of the projectile hadron, and Mfi is the 

effective mass of the fluctuation, which is to be considered on mass shell in the light-cone 
formalism. The kinematical formula for Mji reads, 

,^^m[+fc[^ (4) 

where mj, fcj and aj are the mass, the transverse momentum and the fraction of the total 
light-cone momentum carried by the i-th parton, and we sum over all the partons in the 
fluctuation (X]i«i = !)• 

It is easy to find from (^ the minimal effective mass of a fluctuation containing a cc 
pair, which corresponds to a longest coherence time. 

This is the most general upper restriction for the lifetime of any fluctuation containing a 
cc pair. 

It turns out that at energies corresponding to available data on ^ production we have 
a full coverage of possible scenarios of space-time development. For instance at a;i7' ~ we 
have tc < ^fsjInij^M^. This is only about 1 jvn at E^ = 200 GeV, however, tc becomes 
compatible with the radiuses of heavy nuclei at Tevatron energies. Eh ~ 800 GeV, and it 
is much longer in the latter case ii xp ^ 1. 

It is demonstrated on different exact solutions in 0, ^ ^ that in the approximation 
afJ^{T) ^ 1 the variation of nuclear suppression as function of the coherence length 
always has a form of a linear function of the nuclear longitudinal formfactor squared. 
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where the longitudinal momentum transfer g^ = 1/^c- This formfactor varies from zero in 
the low-energy limit tc <^ Ra up to one at tc ^ Ra- The formula for nuclear suppression. 



which interpolates between the low- and high-energy limits, (|T|) and (0), reads, 

SL = 1 - ^^r (T) [l - Fj(gc)] - ^free (T) FJ(ge) (7) 

In the two following sections we elaborate more with the cross sections CTfree and cfj^. 

3 The freeing cross section and the soft mechanism of charmo- 
nium production 

The soft spectator mechanism of heavy flavour production was the subject of my talk in 
Hirschegg in 1995 [§]. Let me summarize the observations. 

First of all, one should discriminate between the interaction cross section of a hadronic 
fluctuation, containing a heavy quark pair and the cross section of freeing (production 
of) this pair. The former is large, since the fluctuation as a whole is big, while the 
latter is usually small, of the order of the inverse heavy quark mass squared. Indeed, 
soft (low kx) long-wave (in transverse direction) gluons, providing a large cross section, 
cannot resolve such a small-size structure as a heavy-quark fluctuation. The freeing 
amplitude can be represented as a difference between interaction amplitudes of two Fock 
states, which are identical, except one contains the cc fluctuation, but another one does 
not. This is the general principle originating from the pioneering ideas of Feinberg- 
Pomeranchuk and Good- Walker for diffractive production. If all the fluctuations of the 
projectile hadron have the same interaction amplitudes, nothing new can be produced, 
but the initial hadron, which coherence is not disturbed in this case. 

Let us consider a concrete example of a fluctuation of a light quark, which consists of 
the same quark and a colorless cc pair. The freeing amplitude is equal to the difference 
between the amplitude of inelastic interaction of this fluctuation and that of the single 
quark. Although each of these amplitudes is infrared divergent, their difference is finite. 
Since only the light quark-spectator can interact (interaction of the colorless cc pair would 
lead either to open flavour production, or is a higher order l/wig effect) the difference of 
the inelastic amplitudes comes from the difference of the impact parameters of the quark 



with and without the cc fluctuation. If the transverse separation between the quark and 
the center of mass of the cc pair is r-r, the center of gravity of the whole fluctuation (which 
coincides with the impact parameter of the parent hght quark) is distant by (1 — a)rT 
from the cc pair and by ar-p from the hght quark, where a is the fraction of the hght-cone 
momentum of the projectile quark. Thus, the freeing cross section equals to the total 
interaction cross section of a light qq pair with separation ar^. Therefore, at small r-r 
color transparency leads to (Jfreei'^^T) oc r^. This was an important observation of ( see 
also 0). 

The mean transverse separation (tt) between the cc and the q can be estimated as 
follows . The wave function of the fluctuation in the momentum representation is given 
by the energy denominator, ^ccqikr) oc l/{Mji — m^) oc l/{k^ + M^^l — a) + a'^m'^), 
where (§) is used for M?;. Switching to the impact-parameter representation (only for 
the transverse momentum) one gets the wave function of the qcc fluctuation (the spin 
structure is neglected) "^ccqirTyd) oc -K'o(t^t) 0, where r^ = (1 — a)M| + a^m^. This 
distribution gives the mean separation (r^) ~ 1/t^, and the freeing cross section reads, 

„2 
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^free{a) ~ .^ „ ^,, ,^2 , .o 



As different from naive expectations, eq. (^ says that the cross section of freeing of a 
heavy quark fluctuation is small unless the factor (1 — a) is small ~ rri^/M^, i.e. the 
cc takes the main fraction of the initial quark momentum. In this case the freeing cross 
section is as large as a typical hadronic one, ~ 1/m^ (actually, rrig should be treated as 
an effective infra-red cut off for rj- of the order of the inverse conflnement radius, rather 
than the quark mass). This fact may have relevance to the well known problem of a 
substantial increase of nuclear suppression of J/^ production at large xi = p^/p'^, which 
corresponds to an effective absorption cross section close to that for light hadrons. We 
come back to this problem below. 

The production charmonium cross section is given by a convolution of the cross section 
(P) with the momentum distribution function Fh{y) of the quark in the projectile hadron 
(more generally, the distribution of those partons, which participate in the fluctuation 
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containing the cc). 

1 

{(^'freei^l)) '^ J ^^^T J dy {"^ ccg{rT, Xl/y)\'^ Fh{y) (Tfreeirr) 



(9) 
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Assuming the end-point behaviour for Fh(y) oc (1 — y)'^, one arrives to the following 
Xi-dependence of the freeing cross section at large xi, 



W'freeiXl)) OC 
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{in (a, [(1 - xi)M|/m^] +c)} 



(10) 



The double-logarithmic dependence on Ml corresponds to the well known violation of the 
Bjorken scaling in deep-inelastic scattering given by the evolution equations. However, 
the hard scale in (|T0|) is imposed by (1 — xi)M|, rather than M|, prescribed by the 
conventional factorization theorem. Another violation of the factorization is due to the 
constant c in ([T0|), which depends on the exponent (3 in the hadronic structure function. 
Amazingly, the contribution of the soft asymmetric fluctuations with 1 — a ~ my Ml, 
corresponding to large freeing cross section (see (P)), does not vanish at high M|. More- 
over, these soft fluctuations dominate in nuclear shadowing, which is given by the same 
expression (0), but with higher powers of (yfj-f.^.. This can be explained in a simple 
way (like in deep-inelastic scattering [0)if to divide all the fluctuations to two classes, 
the hard one with a free ~ V^l; ^^^ the soft one with a free ~ V'^g- The results for 
the production cress section and the first order shadowing correction are summarized in 
Table 1. 



Table 1. Contributions of soft and hard fluctuations to the charmonium 
production cross section and to nuclear shadowing 



Fluctuation 


\^c-eq ^ 


hN 
"tot 


l^ccq <J tot 


^ccq\ \<Jtot ) 


Hard 


~ 1 


~ 1/A/l 


~1/A/| 


~ 1/M| 


Soft 


~ my Ml 


r^l/ml 


~1/M| 


~ l/mlMl 



One can see from the Table that the soft component dominates the first-order nuclear 
shadowing term proportional to ([o'fjvl'^r)]^), which has the same 1/M| suppression and 



the the impulse approximation term (cr*^(r'p)). This means that shadowing scales in 
M^. This is the common feature of deep- inelastic scattering p], Drell-Yan pair and heavy 
flavour production p , where the soft and hard contributions are known to be of the same 
order. 

Although the dynamics of charmonium production in hadronic collisions is much more 
complicate, one may try to develop a phenomeno logical approach relying on the general 
ideas on space-time development (section 2) and on the dynamics of nuclear attenuation 
(this section). 

Let us classify different contributions to nuclear suppression factor (|^). 

1. Since the dominant contribution to nuclear shadowing comes from asymmetric 
fluctuations q —>■ ccq with a —>■ 1, the J/\l/ has the same xi-distribution as the valence 
quark, i.e. ~ (1 — xi)^ if the projectile is a proton. This includes also the so called color- 
octet mechanism which corresponds to the fluctuation q -^ qcc2g. The exponent in the 
xi-distribution may be even smaller if the colorless cc fluctuation is produces collectively 
by two or three valence quarks, as was suggested in 0. 

Note that in the dominant asymmetric fluctuation 1 — a ~ m'^/M^ and the light 
quark contributes to the effective mass of the fluctuation (|^) the same amount as the cc 
pair. Therefore, the formfactor -Fl(gc) suppressing the soft fluctuations in (0) should be 
evaluated at qc ~ 2g,j,, where g^ = M^/2E\^ is according to (|^) the minimal longitudinal 
momentum transfer. 

The soft freeing cross section as can be borrowed from the data on nuclear shadowing 
in deep-inelastic scattering. In accordance with the modified factorization relation ([ToD we 
use the result of analysis ||T0| of the data on deep- inelastic scattering at Q^ = {1 — xi)M| 
and fix, as ~ 12 mb , assuming xi ~ 0.5. However, at xi > 0.9 (where no data are still 
available) the corresponding value of Q^ is so small that on gets into the domain of the 
vector dominance model, and cr^^ee should be about 20 — 30 mb 

2. Hard mechanisms of charmonium production include direct interaction of the color- 
octet cc fluctuation with the target (e.g. the color-singlet mechanism). Due to color 



screening the freeing cross section is oc (r^g) ^ A/ml ~ 0-0''' f^ (^^^ factor 4 is due to 
color screening suppressing the fluctuations of small size |jll|). Estimated perturbatively 
an ~ 2 m6, substantially smaller than one expects for 0"^*^. 

Since the effective mass of the hard fluctuation is ~ M,j,, the formfactor F\{qc) in (|^) 
should be taken at q^ ~ q^ 

The produced charmonia have the same oc (1 — xiY distribution as the parent gluon. 
Thus, the hard contribution to the freeing cross section is suppressed by (1— xi)" compared 
with the soft mechanism, where n > 2. 

3. The inelastic J/\l/-nucleon cross section is known from experiment (see in [|l[]) with 
a large uncertainty and ranges from 2 to 7 mh. The lowest-order, two gluon graph provides 
af^ ~ 5 ^ 6 mh. 

We expect quite a steep energy dependence of crfj^. Indeed, J /'^ has a smaller radius 
than the light hadrons, and perturbative QCD predicts a growth of the effective Pomeron 
intercept with decreasing radius. The recent measurements of energy dependence of the 
cross section of elastic photoproduction of J/\E' at HERA [^ found that it grows oc s^^^. 



Since this cross section is oc cr*^ ~ [{afJ^Y /\&-kB, one can find the energy dependence of 
a^^^ oc s°'^^, taking into account the energy dependence of the slope parameter B. Note 
that such a high effective Pomeron intercept perfectly agrees with the measured at HERA 
x-dependence of the proton structure function -^^(x, Q"^) at the corresponding virtuality 

The steep energy dependence of af^^ results in a substantial variation of nuclear 
suppression, provided by the first term in eq. (0). 

We expect the same energy-dependence for the hard component, while the soft con- 
tribution is supposed to grow slowly, oc s^'^ . 

4 Evolution of a cc wave packet in a medium. The ^ /^' puzzle 

The low-energy limit ([l|) was written assuming that the charmonium is produced momen- 
tarily at the point of interaction with a bound nucleon. However, a cc pair is produced. 
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rather than the charmonium, and it takes some time to become the final charmonium, 
even if tc -C Ra- The evolution of the cc wave packet is controlled by a different parameter, 
called formation time/length, 

t > ^^^ , (11) 

which is about five times longer than l^. 

The evolution can be solved using either the path integral technique in quark repre- 
sentation [Q or, what is equivalent, in hadronic representation using the coupled-channel 
approach |]12|. In the latter case the two channel approximation, J/\l' and \E'', turns out 
to have a pretty good, about 10 % accuracy. In this case the effective cross section of 
absorption of J/\l/ in (|I|) and (|^ should be replaced by 

^r =^ ^r [l + eRFliqj)] , (12) 

where qf = 1/lf, and F^{q) is defined in (^. Other parameters are defined in |T^, 



R^ ^ 0.25 is the experimentally known ratio of \1/' to J/\l/ production rates, e = — w2/3 
is the ratio of the off-diagonal to the diagonal diffractive "if — N amplitudes evaluated in 



m 



With these parameters one concludes that if // ^ Ra the effective absorption cross 
section (|l^) is only about 0.6 of that for J/\&. 

Another interesting observation of ||12| concerns nuclear suppression of \I''. Naively, 
using the poor man's formula (|I|) one would expect much stronger suppression of the \1/' 
than J/\I'. However, the corresponding effective absorption cross section turns out to be 
quite different from al 



in 5 






^ V2, 



l + ^^i(./) 



(13) 



where r = crf^^ /af^ ~ 7/3 as evaluated in [T^. It is interesting that with these pa- 
rameters the effective absorption cross sections ([T^) and (|1^) are approximately equal, 
provided that // ^ Ra-, i.e. F\{qf) = 1. This observation naturally explains the surpris- 
ing experimental result of nearly the same nuclear suppression of J/\I' and \E''. Note that 
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in the case of the soft mechanism only the hght projectile quark interacts with the target, 
and the structure of the cc pair also does not affect the nuclear attenuation. 

5 Phenomenology of J/^ production off nuclei 

Summarizing the results of previous sections, eq. (^ can be presented in the unitary form, 

^*4 = exp[-ae//(i^,,xi) (T)] , (14) 

where 

ae//(i?,, X,) = iaf/ [l + eR Fl(g^)] [l - FHq^) 

All the notations here were introduced before, except 7, which is the ratio of the weight 
factors of the soft to the hard components. According to Table 1 this factor is suppressed 
by the J/\E' mass, 7 ~ m'^/M^ ^ 0.004, if to treat rriq as a cut off, which is the inverse 
confinement radius ~ Aqcd- Other parameters are also pretty well known, afj^ = 
cr^ (xiEh)'^, where A = 0.25 and the parameter a^ should provide the cross section 
about 5 -T- 6 mb at the energy few tens GeV. Then we expect o"^ ~ 2 mb. The hard cross 
section an = cr'jj (xiEh)'^, as was mentioned, is expected to be smaller than a*^, of the 
order of 2 mb. Then we expect a% ~ 0.7 mb. The exponent n > 2 is not very important, 
we compare with the data at n = 3. 

To elaborate with low-energy data on J/\l' production one should take care of possible 
effects of induced energy loss, which are due to rescattering of the projectile partons in 
the nucleus [jl3|- This is unimportant at the energy Ep = 800 GeV, where the best data 
exist, but may produce an additional suppression at lower energies. We include this effect 



in a rough way through the extra factor K{Eh, xi) to the expression ([1^) 



A-(^,.,0=(if|l)". (16) 

Here Xi = AE/Eh and AE = k (L), where (L) ^ 3Ra/ 4:[1 + F^{qii,)] is the mean free path 
of a parton in the nucleus, where we took into account the Lorentz stretching of the path 
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of a fluctuation. The density of energy loss dE/dz = —n was fixed at the color-triplet 
string tension k = 1 GeV/ fm. The exponent m = 5 in accordance with the measured xi 
dependence of J/\t' production rate in pp collisions. 
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Figure 1: Nuclear suppression of J/\E' production rate as measured in 
iM' lAl '^'^^ calculated with eq. (\l^). See explanations in the text. 



Fitting the data depicted in Fig. 1 on nuclear suppression from the experiments E772 at 



800 GeV PI and NA3 at 200 GeV M] in pA coUisions, as well as the data of the NA38/50 



experiments at 200 {xi ~ 0.25) [|I^ and 450 GeV {xi ~ 0.1) [|T^ we arrived at the following 
best values of the free parameters in (p!5D, 7 = 0.0037 ± 0.0015, a^ = 1.95 ± 0.018 mb 
and 0"^ = 0.64 ± 0.05 mb, which are in a perfect agreement with our pre-evaluations. 
Comparison with the data, in Fig. 1 and Table 2 demonstrates a nice agreement. 



Table 2. Data on pA collisions from the NA50 experiment at 450 GeV |jT^ and from the 
NA38 experiment at 200 GeV [|1^ compared to our results (bottom row) 



450, W 


450, Cu 


450, Al 


450, C 


200, U 


200, W 


200 Cu 


0.67 ±0.08 


0.75 ±0.08 


0.76 ±0.09 


0.85 ±0.1 


0.63 ±0.13 


0.65 ±0.05 


0.75 ±0.17 


0.74 


0.82 


0.87 


0.90 


0.69 


0.71 


0.80 
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Although our results agree with the NA38/50 data within error bars, they are yet 
systematically higher than the nominal experimental values. This is mostly a result of a 
possible inconsistency between the data from the NA3 and NA38/50 experiments. The 
former are higher and have quite small error bars (see Fig. 1). This region of low xi is 
controlled by the values of a^ and a%, which are the free parameters and do not need to 
vary much in order to fit better either of the data. 

6 DifFractive production of charmonia 

Using the data on nuclear shadowing one can estimate the cross section of diffractive 
production of charmonium in hadron-nucleon interaction, i.e. diffractive excitation of 
the projectile hadron to a state containing the charmonium and light hadrons, while the 
target nucleon is intact. 



It is well known since Gribov's paper \TE\ that nuclear shadowing of the total cross 



section is closely related to diffraction. It is a direct consequence of unitarity that the 
forward diffractive cross section can be represented as 



dM' ^^^^ 



dp^^dM^ p^^o 167r ^"^ ^ ^^^^ 

one averages here over the eigen states of interaction, which may be treated as the pro- 
jectile fluctuations. The cross section (|1^ includes also the elastic channel. 

We are interested only in that part of diffraction in (p!?]), which contains a colorless cc 
pair in final state. In this case one should replace the interaction cross sections in ([T7|) 
by the freeing cross sections. Then one arrives at the same expression, which we had 
for nuclear shadowing term in charmonium production (see section 3) when the nuclear 
formfactor saturates {qcRa <^ 1), 

1-S^^= ,^ \^, (T) / dM' ^''^^ 



where da^^/dM'^dp^ is the differential cross section of diffractive production of charmo- 
nium in hN collisions, i.e. the process hN -^ Xi^iN, where the multiparticle state X^ 
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consists of ^ and light hadrons. M is the effective mass of the Xq,, which is restricted by 
the nuclear formfactor at M^ < M^ ~ E^/ Ra- 

As soon as nuclear shadowing of \1' production is known, we can predict the fraction 
of diffraction in the total cross section of \E' production, which reads 

Ml, 
(mAr+Af*)2 

Here B ^ 5 GeV~^ is the slope parameter of the differential diffractive cross section 
hN -^ Xq,N. Its value should be about a half of the slope of elastic NN scattering, since 
only the form factor of the target nucleon contributes to the slope in diffraction to large 
masses. 

The parameter of the parameter as was fixed (using information from deep-inelastic 
scattering) at as = 12 mb in our analysis in the previous section. Since this value well 
describes the data on nuclear suppression at high energy of J/\E', we can safely use it on 
the same footing as the data. 

Thus, we arrive at the estimate, 

6Z^{Mm) ~ 0.12 , (20) 

which is close to the fraction of the diffractive large rapidity gap events observed at HERA 
||19|| . This is not unexpected result, since, as we emphasized, there is much in common 
between deep-inelastic scattering, Drell-Yan reaction and heavy flavor production, espe- 
cially when it concerns the role of soft interactions, which dominate nuclear shadowing 
and diffraction. 

One should not be confused with the fact that the left-hand sides of relations ([T9|) 
and (^) depend on M^, while the right-hand sides are constant. The high-mass tail 
oc 1/M^ of the diffractive M^-distribution leads to a logarithmic gc-dependence of the 
nuclear formfactor, which we neglected and eliminated the formfactor in (|T9|) and (|20|). 



14 



7 Heavy ion collisions 

In this case either the fluctuations, or the formed charmonium propagate through the 
both colliding nucleus. This may sound puzzhng, indeed, we considered above propaga- 
tion of the fluctuations of the projectile proton through the target nucleus. How such a 
fluctuation can propagate through the projectile nucleus, which has the same velocity as 
the projectile nucleon? This is a typical puzzle of the parton model. It is impossible to 
prescribe a fluctuation either to the beam, or to the target in a Lorentz-invariant way. It 
may be done only in a concrete reference frame. 

Thus, nuclear suppression in AB collision is a product of that in pA and pB inter- 
actions. Of course the kinematics should be treated properly, what is easier to do with 
Feynman variable xp = Xi — M^/xiS. Then the nuclear suppression in AB collision reads 

SlBi^p) = S'^j^ixp) S^Bi-xp) (21) 

Using expression ( p!^ - (p!5| ) it easy to calculate the nuclear suppression factor S^^^^xf) 
for J/^, or applying eq. (|13|), for "$'. It is interesting that due to inverse kinematics 
and the strong Xi^-dependence of Spj^{xp) it turns out that Sj^i^xp) < Sj^i^xp), in 
accordance with experimental observation. However, the observed relative suppression of 
"$' is even stronger than given by (|2l|) . This is possible due to final state interaction with 
the produced hadrons (coomovers), which happens at long times when "$' is formed and 
interacts with a large cross section. 

The same coomovers are a plausible source of the additional nuclear suppression of 
j/xj/ production rate in y4i?-collisions, which was found to be especially strong in lead- lead. 
This is too complicate and uncertain problem to be discussed in this note. 

Summarising, we developed for the first time a phenomenological approach to the 
problem of nuclear suppression of charmonium production, based on the same dynamics, 
which is responsible for nuclear shadowing in deep-inelastic scattering and Drell-Yan re- 
action. The key point is an admixture of soft mechanism of heavy flavour production, 
which scales in M^^, dominates at high xp and naturally explains the observed strong 
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nuclear suppression. Many other features observed experimentally are understood as well. 
Our numerical estimates are in a good agreement with the fit to available data. 

The J/\E' — A^ interaction cross section steeply grows with energy due to increasing 
radius of the gluon cloud. As a result, the J/\E' production rate is expected to be very 
much suppressed due to interaction with nucleons at very high energies. We expect 
aZ^ ^ 10 mh at RHIC and aZf ~ 26 mh at LHC, if J/^ is at rest in the cm. 
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earlier in collaboration with Omar Benhar and Andrzej Zieminski, and our numerous 
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